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ABSTRACT 

Context. The quasar B0605-085 (OH 010) shows a hint for probable periodical variability in the radio total flux-density light curves. 
Aims. We study the possible periodicity of B0605-085 in the total flux-density, spectra and opacity changes in order to compare it 
with jet kinematics on parsec scales. 

Methods. We have analyzed archival total flux-density variability at ten frequencies (408 MHz, 4.8 GHz, 6.7 GHz, 8 GHz, 10.7 GHz, 
14.5 GHz, 22 GHz, 37 GHz, 90 GHz, and 230 GHz) together with the archival high-resolution very long baseline interferometry data at 
15 GHz from the MOJAVE monitoring campaign. Using the Fourier transform and discrete autocorrelation methods we have searched 
for periods in the total flux-density light curves. In addition, spectral evolution and changes of the opacity have been analyzed. 
Results. We found a period in multi-frequency total flux-density light curves of 7.9 ± 0.5 yrs. Moreover, a quasi-stationary jet 
component CI follows a prominent helical path on a similar time scale of 8 years. We have also found that the average instantaneous 
speeds of the jet components show a clear helical pattern along the jet with a characteristic scale of 3 mas. Taking into account 
average speeds of jet components, this scale corresponds to a time scale of about 7.7 years. Jet precession can explain the helical 
path of the quasi-stationary jet component CI and the periodical modulation of the total flux-density light curves. We have fitted a 
precession model to the trajectory of the jet component CI, with a viewing angle 0o = 2.6° ± 2.2°, aperture angle of the precession 
cone Q. = 23.9° ± 1.9° and fixed precession period (in the observers frame) P = 7.9 yrs. 
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OO The radio source B0605 -085 (OH 010) is a very bright at tic nuclei. However, only average acceleration of jet components 

S ■ 'f nT-fo^^c^"^ ? T'fnn'^^^^^rin^/nc ' nTc^'^' J" "i i have been discussed in these papers. 

O , of 0.872 (Stickel et al. 1993). B0605-085 was studied in X- t. , , , « . ., / y u. fDr„;r,<: nc<:f 

, . , , , ., ^TTA^TT^T^A . TTf,^ . The to t al flux-density Tadio ligW cuTve s ol B0605 " 08 5 trom 

r — ■ ray and optical wavelengths with CHANDRA and HST tele- ,, . .^ ^ ... , . v. j- a . • i <-.u . / 

O"- , ■' , ^ , ° ,-„„., ,. „ „ University of Michigan Radio Astronomical Observatory (e.g., 

scopes by Sambruna et al. (2004) and IS one of a lew quasars .,, ^ , innn\ u u- ^ e a- ■ u-iv c i ^ 

— .^ ■,■■•, , , , • ■ ir ir T T-. 1- Alleret al. 1999) show hints of periodic variability. Several out- 

C_? which lets have been detected in the X-rays. Very Long-Baseline > . j • .u / i • . i n • j- 

^^ , ^ •' /, ,TT,Ts , • , . , bursts appeared in the source at regular intervals. Periodic vari- 

Interferometry (VLBI) observations revealed a complex struc- , .,., ^ .• i .• i • • ^ • .• . • u- u ■ 

_ , ,..,., . , , • , V , . ability of active galactic nuclei is a fascinating topic which is 



X 



ture of the radio let of the quasar, with multiple bends and . n j . j . .u ^ r\ ^ c ^- ^ 

, / I ', ,T- , , • not well understood at the moment. Only a tew active galax- 

curves at scales from parsecs to kiloparsecs. High-resolution . , . , ^ ., , ..../. ,. .. , 

,,-_-, . ^ , , ,^. . ° ,., les show evidence tor possible periodicity in radio light curves 

space VLBI observations show that the inner let extends in the , r. •. • . i onni ah . i '^r\f\'^ t^ n 7 ^ -^nni 

, , ^ , ... , ^ ^ . , , (e.g., Ratten et al. 2001, Aller et al. 2003, Kelly et al. 2003, 

H south-east direction and at ~0.2 mas it turns to the north-east ^. ,, . , ^nn/i im . . i onn/i i/^ ji » i '^nn£ 

d ,cs . ^f^f^A-. X,,. , . ■ ... Ciaramella et al. 2004, Villata et al. 2004, Kadler et al. 2006, 

. . . (Scottetal. 2004). When observed at centimeter wavelengths the ,^ , . on.- '^nniz <->• . i '^nnn \T■^^ . . i 

^ ^ , . . , , . , , . , Kudryavtseva & Pyatunina 2006, Oian et al. 2007, Villata et al. 

outer part ot the let continues to extend south-east with multiple ^nnni x, • . • • u. • j- c 

^ , , ., -^ , . ... , 2009). Various mechanisms might cause periodic appearance ot 

turns and at kilo-parsec scales It goes eastwards with a south-east .u . • j- i .u u u v i ^ c- ^ 

. . - . /^ . ^nnr^^ r^, , . . J. outbursts in radio wavelength, such as helical movement ot let 

bend at ~3 arcseconds (Cooper et al. 2007). The kinematics ot , , ^ • j o i/^ i u inm \m\ 

r^n^nr ncr ■ .• . , T^ .. . /^^/^.N r, components (e.g., Camenzind & Rrockenberger 1992, Villata 

B0605-085 let was studiedby Kellermann et al. (2004) at 2 cm or.-. ■ innn /^ * . i oaa/in u i 

X- . ,«^..-.rT^ ^ Ti . . . X- & Ratten 1999, Ostorero et al. 2004), shock waves propaga- 

as part ot the MOJAVE 2-cm survey. Based on three epochs ot ,. , ^, . , irim\ *• j- • . ui ,.■ / 

,^ . ,,^„^ 1^^^ . ^«/Ms • tion (e.g., Gomez et al. 1997), accretion disc instabilities (e.g., 

observations (1996, 1999 and 2001), two let components were „ °^ , .nm t u c r> i j oaacn • » 

, , , , ' ,. ,- , .. • , , Honma et al. 1991, Lobanov & Roland 2005), let precession 

found at 1.6 and 3.8 mas distance from the radio core with mod- , e,. ,. , , ^nm n u . i '^r^t^£ n •. . i onin\ 

^^A^ ^-^A A^An ^^-^ (e.g., Stirling et al. 2003, Bach et al. 2006, Britzen et al. 2010), 
erate apparent speeds ot 0.10±0.21c and 0.18+0. 02c, respec- °, • u. u • j- .i j u j ui i 

• . 4,r . . . A-A^^^fAir oo^ .• . or they even might be indirectly caused by a secondary black 

lively. The kinematics ot B0605 -085 and Its components accel- u i . .■ j .u • • ui i u i • .u 

•' ^ hole rotating around the primary super-massive black hole in the 

Send ojfprmt requests to: N.A. Kudryavtseva, nkudryav@mpifr- center (e.g., Lehto & Valtonen 1996, Rieger& Mannheim 2000). 

bonn.mpg.de However, investigation of repeating processes in active galaxies 

* Member of the International Max Planck Research School (IMPRS) is extremely complicated due to various factors. The first com- 

for Astronomy and Astrophysics plicity is that a light curve is usually a mixture of flares which 



N.A. Kudryavtseva et al.: A possible jet precession in the periodic quasar B0605-085 



are most likely not caused solely by one effect. Outbursts might 
be caused by multiple shock waves propagation, changes of the 
viewing angle or other conditions in the jet, such as magnetic 
field or electron density. Looking for periodicity is also com- 
plicated with possible quasi-periodic nature of variability (e.g., 
OJ 287, Kidger 2000 and references therein) and long-term na- 
ture of periods, when the observed timescales are of the order of 
tens of years and really long observations are necessary in order 
to detect such timescales. In this paper we aim to study the prob- 
able periodic nature of the radio total flux-density variability of 
the quasar B0605-085 and to check what might be a possible 
reason for it, investigating the movement of the jet components 
with high-resolution VLBI observations together with an analy- 
sis of spectral indexes and opacity of B0605-085. 

The structure of this paper is organized as follows: in Sect. 2 
we describe the total flux-density variability of B0605-085 and 
focus on periodicity analysis, spectral changes and frequency- 
dependent time lags of the flares. In Sect. 3 parsec-scale jet kine- 
matics of B0605-085 is presented. In Sect. 4 we apply a pre- 
cession model to the trajectory of B0605-085 jet component. 
In Sect. 5 we discuss and present the main conclusions and in 
Sect. 6 we list main results of the paper. 



2. Total flux-density light curves 

2.1. Observations 

For investigating the long-term radio variability of B0605-085 
we use University of Michigan Radio Astronomical Observatory 
(UMRAO) monitoring data at 4.8 GHz, 8 GHz and 14.5 GHz 
(Aller et al. 1999) complemented with the archival data from 
408 MHz to 230 GHz. Table [1] lists observed frequency, time 
span of observations, and references. All references in the table 
refer to the published earlier data. We analyze total flux-density 
variability of the quasar at ten frequencies obtained at University 
of Michigan Radio Astronomical Observatory, Metsahovi Radio 
Astronomical Observatory, Haystack observatory, Bologna in- 
terferometer, Algonquin Radio Observatory, and SEST tele- 
scope. The UMRAO data at 4.8 GHz, 8 GHz and 14.5 GHz since 
March 2001 are unpublished, whereas data from all other tele- 
scopes are archival and have been published earlier. 

The historical light curve is shown in Fig.[T]and spans almost 
40 years of observations. The longest light curve at an individ- 
ual frequency covers more than 30 years at 8 GHz (see Table[T]). 
The individual light curves at 4.8, 8, 14.5, 22, and 37 GHz are 
shown in Fig. 121 The total-flux density variability of B0605-085 
shows hints for a periodical pattern. Four maxima have appeared 
at epochs around 1972.5, 1981.5, 1988.3, and 1995.8 with al- 
most similar time separation. The last three peaks have similar 
brightness at centimeter wavelengths, whereas the first one is 
around 1 Jy brighter. The data at 22 GHz and 37 GHz are less 
frequent and cover shorter time ranges of observations, but it is 
still possible to see that the last peak in 1995-1996 is of similar 
shape and reaches the maximum at the same time as at lower 
frequencies. 

The flares at 14.5 GHz show sub-structure with the dip in 
the middle of a flare of about 0.5 Jy, which is ~30% of flare 
amplitude, suggesting more complex structure of the outburst or 
possible double-structure of the flares. Figure[3]shows a close-up 
view of the complex sub-structure of flares at 14.5 GHz. In order 
to check whether the flares might show the double-peak struc- 
ture we have fitted the sub-structure of outbursts at 14.5 GHz 
with the Gaussian functions (dotted lines in Fig. [3]). The sum of 
the Gaussian functions (solid line) follows the light curve fairly 



well and the possible double peaks have similar separation, they 
have time difference of 2.11 years between 1986.43 and 1988.54 
sub-flares and 1.92 years between the 1994.45 and 1996.37 sub- 
flares. 

We searched also for archival data at optical wavelengths. 
However, due to close proximity of a very bright foreground star 
of 9th magnitude, the optical observations of this quasar are dif- 
ficult and there are not enough data for reconstruction of optical 
light curves. 

2.2. Periodicity analysis 

In order to check for possible periodical behavior of flares in 
the B0605-085 total flux-density radio light curves we have ap- 
plied the Discrete Auto-Correlation Function (DACF) method 
(Edelson & Ki-olik 1988, Hufnagel & Bregman 1992) and the 
date-compensated discrete Fourier transform method (Ferraz- 
Mello 1981). The discrete auto-correlation function method per- 
mits to study the level of auto-correlation in unevenly sam- 
pled data sets avoiding interpolation or addition of artificial 
data points. The values are combined in pairs (a,, bj), for each 
< i, j < N, where N is the number of data points. First, the 
unbinned discrete correlation function is calculated for each pair 
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where a, b are the mean values of the data series, and cr^, <Jh 
are the corresponding standard deviations. The discrete correla- 
tion function values (DCF) for each time range Af^ = tj - tj are 
calculated as an average of all UDCF values, which time inter- 
val fall into the range t - At/2 < Af,j < t + At/2, where r is 
the center of the bin. The higher is the value of At the better 
is the accuracy and the worse is time resolution of the correla- 
tion curve. In case of the auto-correlation function, the signal is 
cross-correlated with itself and a - b,a -b, and cr„ = cri,. The 
error of the DACF is calculated as a standard deviation of the 
DCF value from the group of unbinned UDCF values 



^(^) - WTJ (Z[UDCF,v - DCF{T)ff" . 



(2) 



The DACF method yields different numbers of UDCF per bin, 
which can affect the final correlation curve. In order to check the 
results of the DACF we also applied the Fisher z-transformed 
DACF method, which allows to create data bins with equal num- 
ber of pairs (Alexander 1997). 

The Date-Compensated Discrete Fourier Transform 
(DCDFT) method was created in order to avoid the problem 
of finiteness of the Fourier transform and allows us to estimate 
timescales of variability in unevenly sampled data with better 
precision. The method is based on power spectrum estimation 
fitting sinusoids with various trial frequencies to the data set. In 
case of a periodic process formed by several waves, the DCDFT 
method permits to filter the time series with an already known 
frequency and find additional harmonics. 

We searched for periods in the light curves at 4.8 GHz, 
8 GHz, and 14.5 GHz, which span more than 30 years. We 
have removed a trend before applying a periodicity analysis. The 
trend was estimated by fitting a linear regression function into 
the light curve. The auto-correlation method reveals periods of 
7.7+0.2 yr at 14.5 GHz, 7.2+0.1 yr at 8 GHz, and 8.0+0.2yr at 
4.8 GHz with high correlation coefficients of 0.84, 0.77, and 0.83 
respectively. We have estimated the periods by fitting a Gaussian 
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Table 1. List of radio total flux-density observations used in the paper 



Time Range Telescope 



Reference 



408 MHz 


1975-1990 


Bologna 


4.8 GHz 


1982-2006 


Michigan 


6.7 GHz 


1967-1971 


Algonquin 


8 GHz 


1970-2006 


Haystack, Michigan 


10.7 GHz 


1967-1971 


Algonquin 


14.5 GHz 


1970-2006 


Haystack, Michigan 


22 GHz 


1989-2001 


Metsahovi 


37 GHz 


1989-2001 


Metsahovi 


90 GHz 


1986-1994 


SEST 



230 GHz 1987-1993 SEST 



Bondi et al. 1996a 

Aller et al. 1999, Aller et al. 2003 

Medd et al. 1972 

Dent & Kapitzky 1976, Aller et al. 1999, Aller et al. 2003 

Medd et al. 1972, Andrew et al. 1978 

Dent & Kapitzky 1974, Aller et al. 1999, Aller et al. 2003 

Terasranta et al. 1998 

Terasranta et al. 1998 

Steppe et al. 1988, Steppe et al. 1992, 

Tornikoski et al. 1996, Reuter et al. 1997 

Steppe et al. 1988, Tomikoski et al. 1996 



X 
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408 MHz, Bondi etal. 1996a 

4.8 GHz, Michigan 

6.7 GHz, Medd etal. 1972 ^ 

8 GHz, Michigan and Dent&Kapitzky 1974 

10.7 GHz, Medd etal. 1972 

10.7 GHz, Andrew etal. 1978 

14.5 GHz, Michigan and Dent&Kapitzl<y 1974 

22 GHz, Metsahovi 

37 GHz, Metsahovi 
■•— 90 GHz, Steppe 1988, 1992, Reuter 1997, Tornil<oski 1996 
• 230 GHz, Steppe 1988, Tornil<oski 1996 



1970 



1980 



1990 



2000 



Time [yrs] 



Fig. 1. Historical total flux-density light curves of B0605-085 at ten frequencies from 408 MHz to 230 GHz (see text and TablelTJ. 



function into the DACF peak. The error bars are calculated as 
the error bars from the fit. The Fisher z-transformed method for 
calculating DACF provides the same results. An example of cal- 
culated discrete auto-correlation function at 14.5 GHz is shown 
in Fig. ID (Left). The plot shows a smooth coiTelation curve with 
a large amplitude of about 0.8, which is an evidence for a strong 
period. The date-compensated Fourier transform method also 
shows presence of a strong periodicity of about 8 years. This 
method indicates a period of 8.4 years with an amplitude of 
0.6 Jy at 14.5 GHz, 7.9 years with amplitude of 0.6 Jy at 8 GHz, 
and 8.5 years with amplitude of 0.3 Jy at 4.8 GHz. In Fig. |4] 
(Right) we plot a periodogram calculated for the 8 GHz data. 
The dashed line marks a 5 % probability to get a periodogram 
value by chance. The peak corresponding to the detected period 
is clearly distinguishable from the noise level and has a high am- 
plitude of about 0.6 Jy. The periodogram peak reaches 166.01 at 
the frequency 0.1268 yr ', which coiTesponds to the period of 
7.88 years. The probability to get such peak by chance is less 
than 10 '^. The estimated periods are much longer at 4.8 GHz 
and 14.5 GHz than at 8 GHz, which might be due to different 
length of observations at different frequencies. The light curves 



at 4.8 GHz and 14.5 GHz span about 23 years (if we do not take 
into account a few data points in 1970 at 14.5 GHz), whereas 
the light curve at 8 GHz lasts for 35 years, includes 4 periods 
and therefore provides more reliable estimate of the period. The 
periodicity analysis might show longer periods at 4.8 GHz and 
14.5 GHz with the lack of the first peaks in the light curves. 
Comparison of the obtained sinusoidal signal, where the wave- 
length and phase are equal to the measured period, with the ob- 
served light curves is shown in Figs. |5] (14.5 GHz and 8 GHz) 
and|6](4.8 GHz). It is clearly seen from the plots, that these si- 
nusoids follow fairly well the observed flux-density variability 
for all four flares. The main difference is the absence of the last 
peak in 2003-2006 (expected maximum in 2005.9 at 4.8 GHz, 
2002.9 at 8 GHz, and 2004.1 at 14.5 GHz), which is predicted 
by the periodicity analysis, but is not observed. Both methods 
give similar results and one can calculate an average period of 
7.9 + 0.5 years for all three frequencies and for the two methods. 

The four main peaks in the total flux-density light curve in- 
dicate a periodical behavior of about 8 years. If the period pre- 
serves over time, we would expect a powerful outburst in ~ 2004. 
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Fig. 4. Left: Discrete auto-correlation function for 14.5 GHz measurements; Right: Results of the date-compensated discrete Fourier 
transform of the 8 GHz flux-density data. The dashed line marks a 5 % probability to get a peak by chance. 
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Fig. 5. Total flux-density light curves of B0605-085 at 14.5 GHz and 8 GHz. The red line shows sinusoids derived from the 
date-compensated Fourier transform method with periods of 8.4 and 7.9 years respectively. 



However, after year 2000, the flux-density of B0605-085 stayed 
at the same flux level of about 1 .7 Jy at all five frequencies. 

2.3. Frequency-dependent time delays 

In case of variability caused solely by jet precession we can 
expect that the flares will appear simultaneously at all frequen- 
cies. In order to check whether the peaks at diff'erent frequencies 
reach a maximum at the same time, we calculated frequency- 
dependent time delays. The Gaussian functions were fitted to 
flie light curves at 4.8 GHz, 8 GHz, 14.5 GHz, 22 GHz, 37 GHz, 
and 90 GHz as was described in Pyatunina et al. (2006, 2007). 
The frequency-dependent time delays were estimated as the time 
difference between the Gaussian peaks at different frequencies. 
Due to insufficient data, it was possible to calculate frequency- 
dependent time delays only for the outbursts happened in 1988 
and 1995-1996. However, from Fig. [U it is seen that the 1973 
flare appeared almost simultaneously at 8 GHz and 14.5 GHz. 
We were not able to get reliable fits for the frequencies 22 GHz 
and 37 GHz due to sparse observations. 



The parameters of the Gaussians, fitted to the light curves, 
are shown in Table |2] where frequency, amplitude of a flare, 
time of the maximum, width of a flare ©, and time delay are 
listed. We calculated frequency-dependent time delays with re- 
spect to the position of the peak at the highest frequency (22 GHz 
for the 1995-1996 flare and at 90 GHz for the 1988 flare). It 
is clearly seen that the C outburst in 1988 appeared almost si- 
multaneously at all four frequencies with the negligible differ- 
ence in time between individual frequencies of about 0.08+0.05 
years. The 1995-1996 flare D shows similar properties, the light 
curves at all frequencies have reached the maxima almost si- 
multaneously in about 1995.9, except for 4.8 GHz, which has a 
delay of 0.43+0.07 years. On the other hand, the data at 4.8 GHz 
are poorly sampled during the flare rise, which could shift the 
Gaussian peak. We have also calculated frequency-dependent 
time lags using the cross-correlation function. Table[3]shows de- 
lays between various frequencies for C and D flares. The delays 
for the C flare are zero within the eiTor bars. The D flare has 
shown the large time delay of 0.69 + 0.05 years between 4.8 GHz 
and 22.0 GHz, whereas the lags at all other frequencies are much 
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shorter, from 0.08 + 0.02 to 0.19 + 0.03 years. These values are 
consistent with the lags obtained from the Gaussian fitting, es- 
pecially that we obtained the larger time delay at 4.8 GHz and 
lower lags at other frequencies. The 1982 flare is poorly sam- 
pled, but it is seen from the plot in Fig. [1] that it appeared almost 
simultaneously at 8 GHz and at low frequency 408 MHz. The 
1973 flare was observed only at 8 GHz and 14.5 GHz. There is 
not enough data for this flare to fit Gaussian functions, but it is 
seen from Fig.[T]that rising of the flare appear simultaneously at 
two frequencies. Therefore, we can conclude from the Gaussian 
fitting and from the visual analysis of the flares, that the 1973, 
1988, and 1995-1996 outbursts appeared almost simultaneously 
at different frequencies. 

The bright outbursts appearing simultaneously at all frequen- 
cies can be an evidence for a periodic total flux-density vari- 
ability caused by the jet precession. We would expect that if 
the total flux-density variability is solely due to changes in the 
Doppler factor, then the peaks of flares should appear at the same 
time at various frequencies, since the flux-density S j is chang- 
ing like S j - 5 '.5(0, y)''^"', where a is the spectral index (Lind 
& Blandford 1985). 

Table 2. B0605-085: Parameters of outbursts 



Freq. 


Amplitude 


' max 





Time delay 


[GHz] 


[Jy] 


[yrl 


[yrl 


[yrl 


Cflai-e 










4.8 


1.36 ±0.01 


1988.22 ±0.02 


5.61 ±0.02 


-0.02 ±0.06 


8.0 


1.66 ±0.01 


1988.29 ±0.01 


3.94 ±0.01 


0.06 ± 0.05 


14.5 


1.93 ±0.01 


1988.30 ±0.01 


3.32 ± 0.02 


0.08 ± 0.05 


90.0 


2.21 ± 0.25 


1988.22 ±0.04 


1.71 ±0.05 


0.00 ± 0.04 


D flare 










4.8 


1.30 ±0.01 


1996.24 ± 0.02 


4.51 ±0.02 


0.43 ± 0.07 


8.0 


1.71 ±0.02 


1995.84 ±0.01 


3.33 ± 0.02 


0.03 ± 0.07 


14.5 


1.67 ±0.01 


1995.91 ±0.01 


3.16 ±0.01 


0.09 ± 0.07 


22.0 


1.56 ±0.04 


1995.81 ± 0.05 


2.71 ± 0.05 


0.00 ± 0.05 



Table 3. B0605-085: Time delays estimated with the cross- 
correlation method 



Freq. 1 


Freq. 2 


Time delay 


[GHz] 


[GHz] 


[yr] 


C flare 






4.8 


90.0 


0.12 ±0.19 


8.0 


90.0 


0.34 ± 0.08 


14.5 


90.0 


0.16 ±0.24 


90.0 


90.0 


0.00 


D flare 






4.8 


22.0 


0.69 ± 0.05 


8.0 


22.0 


0.19 ±0.03 


14.5 


22.0 


0.08 ± 0.02 


22.0 


22.0 


0.00 



2.4. Spectral evolution 

In order to check whether all outbursts in the total flux-density 
light curves have similar spectral properties, we constructed 
the quasi-simultaneous spectra for all available frequencies. We 
have selected and plotted spectra for various states of variability 
of B0605-085, such as minimum, rising, maximum and falling 



states. Figure [T] shows the spectral evolution observed for the 
1988 (Left) and 1995-1996 (Right) bright outbursts. The 1988 
flare starts from a steep spectrum in 1983 and 1984, that then 
changes its turnover frequency in 1986. During the maximum 
of the 1988 flare, the spectrum becomes flat and then steepens 
again during the minimum state after the flare. The 1995-1996 
flare shows similar spectral evolution, with a steep spectrum at 
the beginning and at the end of the flare and a flat spectral shape 
during the maximum (Fig.|7]i. 

The light curve of B0605-085 consists of four periodical 
flares in 1973, 1981, 1988, and 1995-1996. The last two flares 
have similar spectral evolution and similar brightness. The ob- 
servations during the 1973 flare are much sparser, but historical 
radio light curve shows that the fluxes at 8 GHz and 14.5 GHz 
are similar, which indicates a flat spectrum. However, it is worth 
to notice that the second 1981 flare shows a completely different 
spectral behavior. One can see from the historical light curve (see 
Fig. [TJ that the spectrum of this flare is steep, since the bright- 
ness at 408 MHz significantly exceeds the brightness at centime- 
ter wavelengths. Therefore, we can conclude that although the 
periodical flares look almost similar in brightness, they reveal 
different spectral properties. 



3. Kinematics of B0605-085 

3.1. Observations 

For an investigation of B0605-085 jet kinemat- 
ics we analyze 9 epochs of VLBA observations at 
15 GHz from the MOJAVE & 2cm survey programmes 
(http://www.physics.purdue.edu/astro/MOJAVE/). The obser- 
vations were performed between 1995.6 and 2005.7. The data 
were fringe-fitted and calibrated before by the MOJAVE collab- 
oration (Lister & Homan 2005, Kellermann et al. 2004). Lister 
et al. (2009) and Homan et al. (2009) discussed jet components 
movement and acceleration in this source in a framework of 
statistical studies of the sample of 135 active galaxies. In this 
paper we aim to study motion of B0605-085 jet components in 
more detail, avoiding pre-established initial models during the 
model-fitting process. 

We performed map cleaning and model fitting of circu- 
lar Gaussian components using the Difinap package (Shepherd 
1997). We used uniform weighting scheme for map cleaning. 
Circular components were chosen in order to avoid extremely 
elongated components and to ease comparison of features and 
their identification. In order to find the optimum set of com- 
ponents and parameters, we fit every data set starting from a 
point-like model, adding a Gaussian component after each run of 
model fitting. The position errors of core separation and position 
angle were estimated as Ar - {dcTrms V^ + 'Speak/o'rms)/25peak 
and A0 = arctan(Ar/r), where (Ttms is the residual noise of the 
map after the subtraction of the model, d is the full width at half 
maximum (FWHM) of the component and 5 peak is the peak flux 
density (Fomalont 1999). However, this formula tends to under- 
estimate the eiTor if the peak flux density is very high or the 
width of a component is small. Therefore, uncertainties have also 
been estimated by comparison of different model fits (±1 com- 
ponent) obtained for the same set of data. This permits to check 
for possible parameter variations of individual jet components 
during model-fitting. The final position error was calculated as a 
maximum value of eiTor bars obtained by two methods. 

The final hybrid maps with fitted Gaussian components are 
shown in Fig. [8] We list all parameters of jet components, such 
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Fig. 7. Spectral evolution of B0605-085 during the 1988 flare (Left) and 1995-1996 (Right). It is clearly seen how spectral shape 
changes from steep at the beginning and at the end of the flares to flat during the maxima. 



as the core separation, the position angle, the flux-density, the 
size, and identification for all epochs in Table |5] 



3.2. The component identification 

The component identification was performed in such way that 
the core separation, the position angle and the flux of jet features 
change smoothly between adjacent epochs. In order to obtain the 
best identification, the jet components were labelled in different 
ways and the identification which gave the smoothest changes 
in all parameters of jet features, such as mean flux and position, 
was chosen as final. The core separations and the flux densities 
of each individual jet component are shown in Fig. |9] left and 
right correspondingly. Individual jet features are marked by var- 
ious colors. Jet feature CI appears quasi-stationary at an average 
core separation of 3.7 mas with a flux density of about 0.1 Jy at 
all epochs. It is marked with a dotted line on the plot. This quasi- 
stationary feature is well-detected in the first four epochs, but af- 
ter 2002 it is blended with the outwards moving jet component 
d2, which was ejected in ~ 1998. This feature may be explained 
with a co-existence of a stationary component and jet features 
moving outwards. Four outwards moving jet components (da, 
dO, dl, d2, and c/3) were ejected during the time of the observa- 
tions, whereas the components d6 and dl are located at a range 
of core separations between 6 and 8 mas, showing slow outward 
movement and have probably been ejected earlier. The jet com- 
ponents which were ejected at the time of observations are at 
a level of 0.3 - 0.4 Jy after ejection and then their flux-density 
fades down to 0. 1 Jy level or less. 



3.3. Component ejections 

We found seven jet components moving outwards with apparent 
speeds from 0. 12 to 0.66 mas/yr, which coiTespond to speeds of 
5.6 c - 3 1 c (see Table|4]i. The d3 component has been ejected in 
1996.9 during the 1995-1996 outburst, after the flare has reached 
its peak. All other components ejected after d3 feature, like dl, 
dO, and da show decreasing speeds for almost each successive 
component. The jet component d3 has a speed of 0.66 + 0.01 
mas/yr, the next feature d2 has speed of 0.57 + 0.01 mas/yr, 
dl has 0.44 + 0.01 mas/yr, dO has 0.46 + 0.02 mas/yr, and the 
last component da has 0.12 + 0.03 mas/yr. A few faint Gaussian 



components could not be identified between different epochs 
(marked by black dots in Fig.|9]l- These components reveal flux- 
densities on the order of 0.05 Jy, and can either be blobs which 
were ejected before 1995 and can not be traced because of lack 
of data, or can be possibly trailing components (e.g., Agudo et 
al. 2001, Kadleretal. 2008). 

The Doppler factor of a jet depends on the Lorentz factor y, 
the angle between the jet flow direction and the line of sight cp, 
and the speed of a jet component in the source frame /3: 



S^ 



I 



y(l -/3cos(p) 



(3) 



Using velocity of the jet component with the highest speed one 
can calculate the maximal viewing angle 



sin0max = 2/?app/(l+/?app) 

and a lower limit for the Lorentz factor 



7a 



a/^ 



app' 



(4) 



(5) 



assuming that the highest apparent speed detected defines the 
lowest possible Lorentz factor of the jet (e.g., Pearson & Zensus 
1986). The fastestjet component of B0605 -085 is d3, which has 
a speed of 31c. The corresponding Lorentz factor is j^in = 31 
and the upper limit for the viewing angle is 0n,ax = 3.7°, which 
gives us an estimation for the Doppler factor 6mm - 12.4. For 
smaller viewing angles (0 -^ 0) the Doppler factors approach a 
value of 62, which yields a range of the Doppler factors 6 = 12.4 
to 62 at viewing angles between 3.7° and 0°. 

The instantaneous speeds of the jet components have been 
calculated as the core separation difference between two subse- 
quent positions of a jet component divided by the time elapsed 
between these two observations. The jet components dl, d3, and 
CI show significant changes of instantaneous speeds which in- 
dicates that they experience acceleration and deceleration along 
the jet (see Fig. [10] for the jet component dl). Figure [TT] shows 
the averaged instantaneous speeds of all jet components in the 
1-mas core separation bins. We have calculated instantaneous 
speeds between each data point for each jet component and then 
have averaged all the instantaneous speeds for data points lo- 
cated within a certain bin of core separation. The average instan- 
taneous speeds of all jet components accelerate at about 0. 1 mas. 
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Fig. 9. Left: Core separation as function of time at 15 GHz. Individual jet components are marked by various colors. Right: Time 
variability of jet components' flux density. 



Table 4. Summary of features' speeds based on the model-fits. 
Component identification, mean core separation, proper motion 
of the component, apparent speed of the component, time of the 
back-extrapolated component ejection, and the maximum view- 
ing angle are listed. 



ID 


'"mean 


/^i- 


/Japp 


to 


0max 




[mas] 


[mas/yr] 


[c] 


[yr] 


[deg] 


dl 


9.26 ± 0.44 


0.20 ±0.10 


9.4 ± 4.7 




12.1 


d6 


7.35 ±0.32 


0.26 ± 0.06 


12.2 ±2.8 




9.4 


d3 


4.31 ±0.67 


0.66 ± 0.01 


31.1 ±0.5 


1996.9 ±0.3 


3.7 


CI 


3.71 ±0.17 


0.09 ± 0.01 


4.2 ± 0.5 




26.8 


d2 


3.11 ±0.56 


0.57 ±0.01 


26.8 ± 0.5 


1997.9 ±0.3 


4.3 


dl 


1.78 ±0.34 


0.44 ± 0.01 


20.7 ± 0.5 


2000.2 ± 0.3 


5.5 


do 


1.33 ±0.22 


0.46 ± 0.02 


21.7 ±0.9 


2001.8 ±0.2 


5.3 


da 


0.39 ± 0.03 


0.12 ±0.03 


5.6 ±1.4 


2001.9 ±0.3 


20.2 



3 mas, 6 mas, and 9 mas distance from the core and decelerate at 
2 mas, 5 mas, and 7 mas, which indicates a characteristic scale 
of 3 mas for acceleration/deceleration of the jet components. The 
helical pattern of the speeds might be due to a helical pattern of 
the jet and multiple bends on the way of the jet components. It 
is worth to notice, that the characteristic spatial scale of 3 mas 
corresponds to a timescale of about 7.7 years (using the average 
speed of the jet components 0.39 mas/yr), which is similar to the 
periodic timescale of variability in the total flux-density. 

3.4. Quasi-Stationary component C\ 

The quasi-stationary component CI is a well-defined feature 
from 1995 to 2001 (see Fig. |9]l- However, after 2002 it is most 
probably blended with the moving jet component d2. Earlier 
than 2002, when the quasi-stationary component is clearly iden- 
tified, there were no moving jet components. That might explain 
why the component CI was blended only after 2002. As the 
newly ejected jet feature approach the position of component CI, 
it becomes more complicated to distinguish between them and it 
is more likely that we see blending of two components. Blending 
and dragging of stationary jet components by moving compo- 
nents was predicted from the magneto-hydrodynamical simula- 



tions by Gomez et al. (1997) and Aloy et al. (2003) and was 
observed in a large number of sources (e.g., 3C 279 Wehrle et al. 
(2001), S5 1803+784 Britzen et al. (2005), 0735+178 Gabuzda 
etal. (1994)). 

The quasi-stationarity of jet component CI can be due to 
various physical mechanisms. It can be caused by jet bending 
or by stationary conical shock waves appearing due to interac- 
tion with the ambient medium (e.g., Alberdi 2000, Gomez et al. 
1997). If the stationary feature in the jet is caused by jet bend- 
ing towards us, then we can expect a significant flux-density in- 
crease in other jet components in the area close to the position of 
the stationary feature due to relativistic boosting. When we plot 
the flux-density versus core separation for all detected Gaussian 
components (excluding the core) in Fig. [12] we see that at core 
separations between 2.5 and 6.2 mas there is a significant in- 
crease in flux-density values. This increase is independent of the 
component identification and might be explained by a jet bend- 
ing in the area around 4 mas. Moreover, this jet bending can ex- 
plain the existence of the quasi-stationary jet feature CI, which 
is located at similar core separation of about 3.7 mas. 

High-resolution VLBI maps of the quasar B0605-085 were 
discussed in several papers (e.g., Padrielli et al. (1986), (1991), 
Bondi et al. (1996b)). A stationary jet component observed at 2 
and 5 GHz in 1997 at a core separation similar to CI was re- 
ported before by Fey&Charlot (2000). They found the jet com- 
ponents at a core separation of 3.0 mas and PA - 100°, which is 
similar to the position of the CI component. This is an additional 
evidence for stationarity of the jet component CI. 

The quasi-stationary jet component CI is located at an av- 
erage core separation of 3.7 mas. However, the position of this 
component changes with time from 2.95 to 4.32 mas. The posi- 
tion angle and the flux-density of CI are also variable. Figure[T3] 
shows these changes over time in core separation, position an- 
gle and flux-density. The position angle and the flux-density of 
the component are anti-correlated. The position angle reaches 
its maximum value of 125 degrees in 2001, when the flux is at 
the minimum value of 0.013+0.002 Jy. There is also a possible 
connection between the core separation and the flux-density of 
the component CI - both of them have minima in about 2001- 
2003 and then experience a rise with a maximum in 2005. In 
rectangular coordinates, the trajectory of the quasi-stationary jet 



N.A. Kudryavtseva et al.: A possible jet precession in the periodic quasar B0605-085 



component CI follows a clear helical path. Figure [14] shows a 
trajectory of CI in the plane of the sky. It is clearly seen from 
the plot that the quasi-stationary component follows a spiral tra- 
jectory. Moreover, the characteristic time scale of a helix turn is 
about 8 years. This time scale is very close to the value of the 
period of the total flux-density variability discussed earlier. This 
might be evidence for a possible connection between the heli- 
cal movement of CI and the total flux-density periodicity. The 
jet might precess and therefore due to variability of the Doppler 
boosting we will see periodic flares in the light curves as well as 
changes of position of a bend in the jet. 

3.5. Comparison of trajectories 

We plot the position angles of all features da-dl, including the 
quasi-stationary component CI versus the core separation in 
Fig. [15] The trajectories of the components d\, d2, d6, dl, and 
CI follow significantly curved trajectories, whereas the trajec- 
tories of the components dO and d3 follow a linear path. This 
can be an evidence for a possible co-existence of two types of 
jet component trajectories, similar to the source 3C 345 (Klare 
et al. 2005). The spread of position angles of all components is 
changing along the jet. In the inner part, within the first 1 mas, 
the position angles are in the range of 123-150 degrees, in the 
middle part of the jet between 1-5 mas, the position angles 
range from 110° to 125°, whereas in the outer part from 5 mas 
to 10 mas, from 95° to 115°. 



since oj ~ l/t. Therefore, we can fit the trajectory in the rectan- 
gular coordinate X(t) of the quasi-stationary jet component CI, 
using the formula [6] We assume that the motion of CI reflects 
the movement of the jet. The non-linear least squares method 
was used for the precession model fitting. The angular velocity 
of the precession w and the redshift of the source z are known 
and we fix them during the fit, whereas the aperture angle of the 
precessing cone D., (po, and t/q are used as free parameters. The 
results of the fit are shown in Fig.[T7](Left) and the parameters of 
the precession model for the parsec jet of B0605-085 are listed 
in Table [6] The precession model fits well the trajectory of the 
quasi-stationary jet component CI. The same parameters were 
used to describe the motion of CI in declination. The model for 
the declination is shown in Fig. [17] (Right). 

We found another set of parameters Q = 92.2°, 0o = 69.8°, 
7/0 = 12.9° as one of the possible solutions. However, the view- 
ing angle is very large. Using the equation 



A 



y6sin( 



app 



1 -yScOSf 



(9) 



16.5c, 



and taking average value of apparent speeds /Sapp.aver 
we can calculate the speed in the source frame /? = 2.49c, which 
is higher then the speed of light. Precession cone aperture angle 
Q. = 92.2° is also unlikely. Therefore, this set of parameters is 
not plausible. The solution in Table|6]is the only physically fea- 
sible, since all the parameter space was covered during the fitting 
procedure. 



4. Precession model 

The probable periodical variability of the radio total flux- 
densities of B0605-085 and the helical motion of the quasi- 
stationary jet component CI can be explained by jet precession. 
We used a precession model, described in Abraham & CaiTara 
(1998) and Caproni & Abraham (2004) for fitting the helical path 
of the jet component CI. In this model, the rectangular coordi- 
nates of a precession cone in the source frame are changing with 
time f due to precession: 



Table 6. Precession model parameters (see Sect.lUi. 



Parameter Description 



Value 



p 


Precession period, fixed 


7.9 yr 


Q. 


Aperture angle, variable 


23.9° ± 1.9° 


00 


Viewing angle, variable 


2.6° ± 2.2° 


TlO 


Projection angle, variable 


33.6° ± 6.5° 



X(t') = [cos Q sin 0o + sin Q. cos (po sin ojt'] cos tjq - 

sinQcosojf' sin 770, 



Y(t') = [cos Q sin 0o + sin Q. cos (po sin ojt'] sin 770 + 

sin D. cos LiJt' cos 770, 



(6) 



(7) 



where Q. is the semi-aperture angle of the precession cone, (po 
is the angle between the precession cone axis and the line of 
sight and it is actually the average viewing angle of the source, 
and ?7o is the projected angle of the cone axis onto the plane of 
the sky (see Fig. [T6] l. The angular velocity of the precession is 
(L) = l/P, where P is the period. We take a P value of 7.9-year 
(timescale found in the total flux-density radio light curves and 
in the helical movement of CI). The time in the source frame 
(f) and the frame of the observer (f) are related by the Doppler 
factor 6 as 



At' = --!—^At. 

l+z 



(8) 



However, the only time-dependent term in the equations ((6]l 
and (O is cjt, which does not depend on this time corrections 



5. Discussion 

Assuming that the helical motion of the quasi-stationary jet com- 
ponent CI is caused solely by jet precession (with a period of 8 
years), we were able to determine the geometrical parameters 
of the jet precession, such as the aperture angle of a preces- 
sion cone, the viewing angle, and the projection angle. The fitted 
viewing angle (po = 2.6° + 2.2° agrees well with the upper limit 
for the viewing angle deduced from the apparent speeds of the 
jet components ^^ax - 3.7°. It is also similar to the viewing an- 
gle (^var -5° derived from the radio total flux-density variability 
by Hovatta et al. (2008). The jet precession will also cause the 
periodical variability of the radio total flux-density light curves 
with the same period of 8 years, due to periodical changes of the 
Doppler factor. The flux density S j is changing with the Doppler 
factor due to the jet precession S j = S'M(l),yy^", where a is 
the spectral index (Lind & Blandford 1985). From this equation 
we can expect that the flares due to the jet precession will ap- 
pear simultaneously at different frequencies, which is actually 
observed for the major outbursts in 1973, 1988, and 1995-1996. 
Therefore, the jet precession is probably responsible for the to- 
tal flux-density variability, periodicity in the outbursts, and the 
helical motion of the quasi-stationary jet component CI. 



N.A. Kudryavtseva et al.: A possible jet precession in the periodic quasar B0605-085 



4.6 

„ 4.4 

w 

n! 4.2 

E 

" 4.0 

■2 3.8 

2 3.6 
ro 

" 3.2 

5 3.0 



0605-085, CI component 



0605-085, CI component 




2.6 ' 

1994 19 



0605-085, CI component 




1998 2000 2002 2004 2006 2008 



Time [yrs] 



1994 1996 1998 2000 2002 2004 2006 2008 
Time [yrs] 



0.00 ' ' ' ' ' ' ' — 

1994 1996 1998 2000 2002 2004 2006 2008 



Time [yrs] 



Fig. 13. Time evolution of core separation (Left), position angle (Middle) and flux-density (Right) of the quasi-stationary component 
CI. 
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Fig. 14. Trajectory of the quasi-stationary component CI in rectangular coordinates in three-dimensional space with the third axis 
showing the time (Left) and in two dimensions (Right). The trajectory of the component seems to follow a helical path. 



Some questions remain unclear. It is still not known why the 
four major outbursts in the total flux-density light curve which 
repeat periodically with an 8-year period have dififerent spectral 
properties. Longer monitoring of this quasar at several frequen- 
cies is necessary in order to understand spectral properties of the 
flares. 

The precession of B0605-085 jet can be caused by sev- 
eral physical mechanisms. Jet instabilities, a secondary black 
holes rotating around the primary black hole, and accretion disc 
instabilities might be responsible for the observed changes of 
jet direction (e.g., Camenzind & Krockenberger 1992, Hardee 
& Norman 1988, Lobanov & Roland 2005). At the moment it 
is impossible to claim exactly which mechanism can produce 
the observed jet precession. It is worth to notice, however, that 
the possible double-peak structure of the flares at 14.5 GHz in 
B0605-085 is similar to the double-peak flares of OJ 287 ob- 
served in optical light curves which were explained by the or- 
biting secondary black hole. At the moment we can not say 
whether there is a secondary black hole in the center of a hosting 
galaxy in B0605-085 and further observations and theoretical 
modelling are needed. 

The next total flux-density flare of the source might appear 
in 2012 if the 8-year period will preserve over time. The 8-year 



period in the flux-density variability has predicted the next pow- 
erful outburst to appear in ~ 2004. However, the flare was of 
a very low flux level which might be due to long-term trend 
in flux-density variability. Recent observations of B0605-085 
have shown a rise in fluxes which might suggest a probable peak 
close to 2012 taking into account the duration of the outbursts 
in this source. This will mean that the 8-year period interferes 
with long-term changes of the flux-densities. Further radio total 
flux-density and VLBI observations are needed in order to check 
this prediction. 



6. Summary 

Our analysis of the jet kinematics in B0605-085 shows: 

- We find strong evidence for a period of the order of 8 years 
(7.9 + 0.5 years when averaged over frequencies) in the total 
flux-density Hght curves at 4.8 GHz, 8 GHz and 14.5 GHz, 
which was observed for four cycles. 

- The measurement of frequency-dependent time delays for 
the brightest outbursts in 1988, and 1995-1996 and the visual 
analysis for the 1973 flare shows that these flares appeared 
almost simultaneously at all frequencies. 
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Fig. 17. Sky trajectory of the quasi-stationary component CI in right ascension (Left) and declination (Right). The solid line shows 
the fit of the precession model. The precession model fits well the trajectory of a quasi-stationary jet component CI 



- The spectral analysis of the flares shows that the main flares 
in 1973, 1988, 1995-1996 have a flat spectrum, whereas the 
flare in 1981 has a steep spectrum, what suggests that the 
four cycles of the 8-year period have different spectral prop- 
erties. 

- The average instantaneous speeds of the jet components re- 
veal a helical pattern along the jet axis with a characteris- 
tic scale of 3 mas. This scale corresponds to a time scale of 
about 7.7 years, which is similar to the periodicity timescale 
found in the total flux-densities. 

- The quasi-stationary jet component CI follows a helical path 
with a period of 8 years. This time scale coincides with the 
period found in the total flux-density light curves. 

- We found evidence that the jet components of B0605-085 
follow two different types of trajectories. Some components 
move along the straight lines, whereas other components fol- 
low significantly curved paths. 

- The fit of the precession model to the trajectory of the quasi- 
stationary jet component CI has shown that it can be ex- 
plained by changes of the jet direction with a precession 
period of P = 8 yrs (in the observers frame), aperture an- 
gle of the precession cone Q. = 23.9° + 1.9°, viewing angle 
^0 = 2.6° + 2.2°, and projection angle 770 = 33.6° + 6.5°. 
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Fig. 2. Total radio flux-density light curves of B0605-085 at 
five frequencies. The solid lines mark positions of the brightest 
flares. 
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Fig. 3. Complex sub-structure of flares at 14.5 GHz. The dotted 
lines show the Gaussian functions fitted to the sub-structure of 
flares. The solid line is the sum of fitted Gaussian functions. 
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Fig, 6, Total flux density light curve of B0605-085 at 4.8 GHz. 
The red line shows a sinusoid with a period of 8.5 years, derived 
from the date-compensated Fourier transform method 
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Fig. 8. Contour images of B0605-085 and overlapped on them 
circular Gaussian components to model the observed visibilities 
(see Table |5]). The distance between images on the plot is pro- 
portional to time passed. 
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Fig. 10. Instantaneous speeds of the jet component dl, which 
experiences significant acceleration when moving along the jet. 
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Fig. 11. Instantaneous speeds of all jet components of 
B0605-085 averaged in 1 mas core separation bins. We inter- 
pret that the jet is helical with a characteristic scale of 3 mas. 
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Fig. 12. Flux-density of all Gaussian components found in the 
jet of B0605-085 at 15 GHz. The core was excluded from the 
plot for smaller range of flux densities. It is clearly seen that all 
components become significantly brighter when they are at the 
range of the core separations from 2.5 to 6.2 mas, which can be 
explained by the jet bending. 
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Fig. 15. Position angle changes of the jet components da-dl, 
including the quasi-stationary component CI at 15 GHz. 
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Fig. 16. Jet precession geometry. The line of sight is parallel to 
the z-axis. 
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Table 5. Model-fit results for B0605-085 at 15 GHz. We list the epoch of observation, the jet component identification, the flux-density, the radial distance of the 
component center from the center of the map, the position angle of the center of the component, the FWHM major axis of the circular component, and the position 
angle of the major axis of the component. 



Epoch 


Id. 


S 


r 


e 


Ma.A. 


[yr] 




[Jy] 


[mas] 


[deg] 


[mas] 




1995.57 


D 


2.533±0.382 


0.0±0.04 


0.0 


0.32 






0.045±0.015 


1.35±0.06 


103.59±4.51 


0.81 




C\ 


0.084±0.014 


3.20±0.02 


108.79±0.41 


0.33 


1996.82 


D 


2.001±0.319 


0.0±0.16 


0.0 


0.55 






0.045±0.022 


1.70±0.07 


110.22±2.05 


0.20 






0.045±0.007 


2.82±0.04 


125.70±0.90 


0.30 




C\ 


0.067±0.011 


3.19±0.05 


107.67±0.97 


0.53 


1998.83 


D 


0.770±0.107 


0.0±0.03 


0.0 


0.10 




dl 


0.281±0.035 


0.50±0.01 


123.65±1.92 


0.43 




(13 


0.248±0.047 


1.38±0.05 


127.61±0.86 


0.51 




CI 


0.099±0.014 


3.80±0.04 


108.49±0.86 


0.76 




d6 


0.022±0.009 


6.43±0.37 


112.29±3.20 


2.59 




dl 


0.008±0.001 


8.38±0.43 


98.75±3.01 


0.84 






0.011±0.002 


9.95±0.98 


108.37±5.50 


1.66 


2001.20 


D 


0.978±0.144 


0.0±0.02 


0.0 


0.11 




dl 


0.384±0.008 


0.43±0.01 


126.68±1.09 


0.39 




dl 


0.073±0.028 


1.69±0.07 


130.95±1.01 


0.97 




d3 


0.058±0.044 


2.32±0.08 


117.61±1.18 


1.08 




CI 


0.013±0.002 


3.61±0.17 


126.00±2.70 


0.61 






0.109±0.013 


4.22±0.05 


112.50±0.76 


1.34 




db 


0.011 ±0.002 


6.71±0.27 


117.50±2.30 


0.87 




dl 


0.013±0.002 


8.83±0.29 


104.45±1.90 


1.21 


2003.09 


D 


1.145±0.167 


0.0±0.02 


0.0 


0.22 




dO 


0.048±0.007 


0.60±0.02 


151.77±1.70 


0.21 




dl 


0.150±0.042 


1.24±0.10 


129.60±1.99 


1.01 




dl + Cl 


0.039±0.007 


2.95±0.08 


123.91±1.40 


0.77 




d3 


0.051 ±0.044 


4.25±0.08 


110.50±2.68 


2.14 


2004.44 


D 


0.843±0.119 


0.0±0.07 


0.0 


0.08 




da 


0.376±0.067 


0.30±0.03 


133.42±3.42 


0.15 




dO 


0.141 ±0.022 


1.31±0.06 


133.02±2.01 


1.17 




dl 


0.031 ±0.003 


1.79±0.02 


110.46±0.63 


0.10 




dl + Cl 


0.096±0.014 


3.72±0.12 


118.60±1.11 


1.98 




d3 


0.058±0.002 


5.41±0.09 


113.42±0.73 


1.18 




db 


0.009±0.001 


7.85±0.18 


118.32±1.30 


0.56 




dl 


0.019±0.010 


9.42±0.39 


108.10±6.79 


3.76 


2005.17 


D 


0.835±0.119 


0.0±0.06 


0.0 


0.05 




da 


0.323±0.106 


0.40±0.06 


133.26±5.41 


0.08 




dO 


0.061 ±0.009 


1.20±0.17 


134.21±2.60 


0.74 




dl 


0.081±0.011 


2.10±0.07 


121.01±1.27 


1.01 




dl + Cl 


0.100±0.016 


4.28±0.06 


121.57±0.87 


1.50 




d3 


0.033±0.005 


5.75±0.08 


114.90±0.80 


0.75 






0.007±0.002 


6.18±0.09 


106.30±0.90 


0.13 




d6 


0.027±0.006 


7.73±0.33 


103.00±2.50 


2.94 






0.005±0.001 


11.73±0.34 


103.70±1.70 


0.48 






0.007±0.001 


13.26±0.31 


101. 00± 1.40 


0.49 






0.012±0.002 


19.83±1.24 


111.90±3.60 


3.99 


2005.56 


D 


0.907±0.132 


0.0±0.03 


0.0 


0.09 




da 


0.358±0.058 


0.44±0.02 


135.27±1.18 


0.16 




dO 


0.084±0.017 


1.77±0.05 


126.35±1.53 


1.13 




dl 


0.019±0.006 


2.54±0.08 


116.81±1.83 


0.42 




dl + Cl 


0.103±0.021 


4.32±0.15 


121.36±1.88 


1.38 




d3 


0.058±0.007 


5.75±0.13 


112.87±0.77 


1.17 




db 


0.009±0.001 


8.04±0.55 


104.90±4.00 


1.68 




dl 


0.004±0.001 


10.41±1.19 


95.45±6.50 


1.46 






0.008±0.001 


11.12±0.43 


110.02±2.20 


1.05 






0.003±0.001 


21.24±3.66 


113.30±9.80 


3.38 


2005.71 


D 


0.822±0.123 


0.0±0.04 


0.0 


0.13 




da 


0.288±0.044 


0.44±0.02 


132.92±1.41 


0.22 




dO 


0.086±0.011 


1.77±0.04 


127.87±1.42 


1.01 




dl 


0.026±0.004 


2.61±0.05 


118.97±1.20 


0.43 




dl + Cl 


0.053±0.008 


4.28±0.06 


124.57±1.18 


1.08 




d3 


0.108±0.016 


5.31±0.09 


114.02±0.98 


2.10 



